Theoretical evidence for a reentrant phase diagram in ortho-para mixtures of solid H 2 
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We develop a multi order parameter mean-field formalism for systems of coupled quantum rotors. 
The scheme is developed to account for systems where ortho-para distinction is valid. We apply 
our formalism to solid H2 and D2. We find an anomalous reentrant orientational phase transition 
for both systems at thermal equilibrium. The correlation functions of the order parameter indicate 
short-range order at low temperatures. As temperature is increased the correlation increases along 
the phase boundary. We also find that even extremely small odd- J concentrations (1%) can trigger 
short-range orientational ordering. 
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Quantum effects dominate the low temperature (T < 
200K) phase diagram of solid molecular hydrogen in 
a wide range of pressures from ambient up to ~ 100 
GPa0,|2. In this regime the coupling between molecules 
is smaller than the molecular rotational constant, so 
quantum effects are generally described by means of 
weakly coupled quantum-rotor models [3|. Homonuclear 
molecules (H2 and D2) can assume only even or odd val- 
ues of the rotational quantum number J, depending on 
the parity of the nuclear spin. Important differences exist 
in the phase diagrams of even- J (para-H2 and ortho-D2), 
odd- J (ortho-112 and para-T>2) and all- J (HD) species. 
(For a summary of experimental results on the orienta- 
tional ordering in H 2 , D 2 , and HD see Fig. lb of ref- 
erence^.) At low pressure or high temperature, even- J 
species are found in a rotationally disordered free-rotor 
state (phase I). Increasing pressure causes an increase 
of the intermolecular coupling, and eventually leads to 
an orientationally ordered state (phase II). Odd- J sys- 
tems on the other hand are orientationally ordered at 
low temperature and ambient pressure and remain or- 
dered as pressure is increased. The stronger tendency of 
ortho-H2 to order can be traced to the fact that its J = 1 
lowest rotational state allows for a spherically asymmet- 
ric ground state, unlike the J = ground state of even- J 
species. The pressure-temperature phase diagram of HD 
exhibits a peculiar reentrant shape ,5] . Reentrance refers 
to phase diagrams where in some range of pressure the 
system reenters the disordered phase at ultra-low temper- 
atures (see HD in Fig. 0. The zero-temperature orien- 
tationally disordered phase is characterized by an energy 
gap against J = 1 excitations. When this gap is suffi- 
ciently small and the temperature is finite, the thermally 
generated J = 1 excitations suffice to induce ordering, 
which is then reentrant, as also shown by mean-field thc- 
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FIG. 1: Mean- field phase diagrams for coupled quadrupolar 
rotor models corresponding to solid molecular hydrogen and 
various isotopes. 
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ory |y, |Y|, |2| . Reentrance is also found in models of two- 
dimensional rotors 0, such as the quan tum anisotropic 
planar rotor (QAPR) model E3, Ell Il3 

High pressure experiments have focused mainly on 
the behavior of pure species. While understanding the 
orientational transition in ortho-para mixtures of quan- 
tum rotors in two and three dimensions is a prob- 
lem of l ong -standing theoretical and experimental inter- 
est [H El El El E3, 111 El 13 EU , less attention has 
been given at establishing the effects of ortho-para dis- 
tinction. Exceptions [l3l Eil EH Eif are the investigation 
of vibrons in ortho-para mixtures by Feldman et al. EH 
and a recent work by Goncharov et al. E3 > where ortho- 
D2 mixed with small amounts of para-Y>2 indicated that 
the possibility of an orientationally frustrated phase be- 
tween phases I and II (phase IF). Phase IF persisted for 
a narrow pressure range (~ 2GPa) for a thermally equi- 
librated ortho-para mixture. 

At ambient pressure hydrogen is known to enter gradu- 
ally a phase of short-range order (orientational glass) ES 
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FIG. 2: Multi order parameter (MOP) mean-field phase di- 
agrams for coupled quadrupolar rotor models corresponding 
to solid molecular hydrogen at various ortho concentrations. 



|2l| when the concentration of ortho species is below ~ 53 
% ^3 13 113 ■ Short-range order is believed to be a 
manifestation of the orientational frustration due to the 
random location of ortho molecules in the lattice, and 
eventually transforms into long-range order at ortho con- 
centrations approaching unity. The pressure-induced in- 
crease of the intermolecular anisotropic coupling causes 
a decrease of the minimum concentration of ortho species 
required to trigger the orientational transition jlj . At ex- 
tremely high pressures (~ 110 GPa) orientational freez- 
ing is observed also in pure para-^2, as stated above. 
Ambient pressure experiments on pure species as well as 
on samples with fixed concentration are possible because 
the rate of interconversion between para and ortho H2 
species is extremely slow at low pressure. Samples with 
fixed ortho-para concentration are prepared by letting the 
system thermalize at the corresponding temperature for 
very long times and experiments are then performed at 
any other temperature by rapid heating or quench ( rapid 
on the time scale of interconversion). At high pressures 
however the interconversion rate rises steeply j^. At 
58 GPa, the highest pressures where the conversion rate 
has been measured, interconversion takes place in the 
time scale of seconds. While interconversion is still or- 
ders of magnitude slower than molecular rotation, the 
ortho-para distribution in this range of pressure is likely 
to be determined by thermal equilibrium. On the other 
hand ortho-para distinction is still valid. 

Previous theoretical studies addressing quantum ef- 
fects in compressed hydro gen (mean- field HSU or 
quantum Monte Carlo |23l \24. I25L l2rj| studies) have in 
most cases ne glect ed the distinction between even- J and 
odd- J species [23j, and in those cases where the distinc- 
tion was explicitly considered, the analysis was restricted 
to pure even- J or odd- J species 0> HE Hi]- On the 
other hand, theoretical modeling of mixed ortho-para sys- 



tems 0, 0, |2(| at low pressure is generally limited to 
J = 0/ J = 1 states. [l3,|2(j This is not an accurate model 
for compressed H2, as pressure and temperature cause a 
non-negligible admixture of rotational states with higher 
J's in the ground state of both ortho and para species. 
A binary alloy model 0] was also used to analyze the 
low-pressure vibrons in ortho-para mixtures. 

In order to calculate the phase diagram of ortho-para 
mixtures of hydrogen in the solid phase we develop a 
multi order parameter (MOP) mean-field theory which 
accounts explicitly for ortho-para distinction in systems 
of coupled rotors. Our formalism enables us to treat a 
system of ~ 3000 molecules, which would be extremely 
computationally intensive with quantum Monte Carlo. 
Our formalism also includes correlations due to ortho- 
para distinction, a feature that is absent in standard 
mean-field theory. We apply the formalism to a system of 
coupled quadrupolar rotors whose centers of mass form 
a face-centered cubic lattice. We find that the phase line 
separating the orientationally disordered state (phase- 
I) from the orientationally ordered state is reentrant in 
the case of thermally equilibrated ortho-para mixtures 
for both H2 and D2. We also find that orientational or- 
der is short-ranged at low temperatures. The validity of 
our conclusion is supported by the fact that for the two- 
dimensional model mean-field theory J^j. Ilfl| is in agree- 
ment with quantum Monte Carlo [111 Il2| simulations on 
the overall shape of the phase diagram. 

The Hamiltonian of a system of N coupled quadrupo- 
lar quantum rotors interacting through a quadrupolc- 
quadrupole potential can be written in the form 

i=l i<j v %3 ' 

J2C(22A;mn)Y 2m (n i )Y 2n (n j )Yi m+n (n ij ), 

where B is the molecular rotational constant, K is the 
coupling strength, C (224; ran) are Clebsch-Gordan co- 
efficients, tti denote the coordinates of rotor i, and f2y 
denote the direction of the vector connecting rotors i and 
j. We define Rq to be the nearest neighbor distance, R\j 
is the distance between molecules i and j. Our MOP 
mean-field theory is based on the trial Hamiltonian 

N N , „ v 5 

H Q = B^U ^Etf 1 ) ( 2 ) 

i— 1 i<j ^ l 3 * 

C(224; 00)1-20(^)7^ 

where 7^ are parameters. Variation of the free-energy 
leads to the self-consistent expression 

7* = (*2o(fii))o. (3) 

Since the trial Hamiltonian in Eq. © is a sum of single- 
rotor Hamiltonians, it follows that for each ji in Eq. © 
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FIG. 3: Correlation functions along the phase boundary of 
the 10% ortho system of solid H2. 

the average needs to be performed over the correspond- 
ing coordinate fij only. Ortho-para distinction can be im- 
plemented by restricting a particular average to be over 
odd- J or even- J states. 

Since in this study we are only interested in the phase 
diagram we expand to first-order in the potential [9j. 
Such an expansion is expected to be valid here, since 
the transition between the long-range ordered state and 
the disordered state is only weakly first-order Fur- 
thermore the transition between the short-range ordered 
state and the disordered state is known to be continuous 
at ambient pressure |l7|. The resulting expression can 
be written 

^^E7,(|)V 4 o(^), (4) 

where $i is the phase correlator defined as <&i = 

Jo dr(Y 2 0(Q(0))Y2o(Q(T))}fre e ,i- The averaging is to be 
performed over a free rotor at inverse temperature /3, and 
ortho-para distinction can be invoked by restricting the 
average as described above. Thus, for a given tempera- 
ture and rotational constant, $j can take on two values 
depending on whether rotor i is odd- J or even- J. Order 
will be signalled by non-zero solutions for 7$, disorder by 
solutions in which all 7, are identically zero. The stan- 
dard mean-field theory 0, Q is automatically recovered 
for pure systems. 

We calculated the phase diagrams for systems of fixed 
odd- J fraction, as well as at the thermal equilibrium dis- 
tribution, taking account of nuclear spin degeneracy. D2 
and H2 differ by the values of the rotational constants 
{2Bjj 2 — Bh 2 ) and in the degeneracies of states. Or- 
tho(para) H 2 is restricted to be odd- J(even- J) angular 
momentum and in D2 the reverse. In our calculations 
for mixed systems, configurations of 2916 molecules on a 
face-centered cubic lattice were used. For a given odd- J 
fraction 20 configurations are generated. A cutoff dis- 



FIG. 4: Phase diagrams of Fig. and those of the pure 
systems H2 and D2 at thermal equilibrium distribution calcu- 
lated using our multi order parameter (MOP) mean-field the- 
ory and the standard mean-field theory (MF). For the mean- 
field phase diagrams see legend of Fig. Q 



tance of 2.65 times the lattice constant is used, thus 134 
nearest neighbors are included in the interaction calcu- 
lation. The phase transition point is obtained by diago- 
nalizing Eq. Q (which is a matrix equation) and finding 
the lowest coupling constant for which a nonzero solution 
exists. 

In Fig. ^we present the results of standard mean-field 
theory phase diagrams of pure odd-J/even-J of H2 and 
D2, and HD, all in excellent qualitative agreement with 
the experimental results (see Fig. lb of reference 4). The 
energy scale is defined to be the rotational constant of 
the H2 molecule. The main difference between the phase 
diagrams of odd- J and even- J systems is accurately cap- 
tured, namely that at low pressures odd- J systems are 
always ordered, whereas even- J systems order at finite 
pressures. As expected from experiment odd-J D2 or- 
ders at a lower coupling strength than odd-J H2, and 
the reentrant phase transition in HD is also well repro- 
duced by mean- field theory [E 0> HI • 

In Fig. the results of the formalism presented above 
are shown for solid H2. As the ortho concentration is de- 
creased the system tends towards disorder entering the 
ordered state at higher coupling constants for a given 
temperature. A noteworthy result of our calculations is 
that even at an ortho concentration of 1% the system 
enters an ordered state at coupling constants quite dif- 
ferent from that of pure para hydrogen, and that for any 
ortho concentration the ground state is always ordered. 
The 50% ortho system is very close to the pure ortho one. 
Our results are consistent with the phase diagram shown 
in Fig. 1 of Ref. in which the pure para-H2 phase 
transition occurs at ~ 110 GPa, and ortho containing 
samples show ordering at lower pressures. 

In order to assess the nature of the ordering we calcu- 
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FIG. 5: Correlation functions along the phase boundary of 
the of solid H2 at thermal equlibrium. 

lated correlation functions of the local order parameters 
{7^}. This is an advantage of our multi order parame- 
ter formalism over the standard mean-field theory, since 
in single-site mean-field theory only phases of complete 
order or disorder are possible. In particular we have cal- 
culated 

G(r) = <( 7 (0)7(0) 2 >- (5) 

The results along the phase line of the 10% ortho sys- 
tem of hydrogen are shown in Fig. |3| The correlation 
increases along the phase boundary with increasing tem- 
perature indicating the onset of long-range order. The 
onset of long range order is due to the fact that ortho-para 
distinction ceases as temperature and coupling constant 
(pressure) are increased. 

The orientational ordering of the system with ortho- 
para concentrations at thermal equilibrium is shown in 
Fig. 01 along with the phase diagrams of the pure systems 
(also shown in Fig. ^| for comparison. For H2 and D2 we 
calculated the phase diagrams via the multi order param- 
eter mean-field theory presented here, and via the stan- 
dard mean-field theory. The resulting phase diagrams 
indicate reentrance in both D2 and H2. Reentrance is 
also seen [i3| m the corresponding 2D model 0. Reen- 
trance is stronger in the case of H 2 , due to higher rela- 
tive weight of odd- J contribution (nuclear spin degener- 
acy), and higher rotational constant (stronger quantum 
effects). The correlation functions for different tempera- 
tures along the reentrant phase diagram are shown for H2 
in Fig. As the temperature increases correlation in- 
creases along the phase boundary. At high temperatures 
(T > 1) the order is definitely long-range. We conclude 
that short-range order may be present up to T = 0.75. 
In the case of D2 reentrance is less severe. 

In summary we calculated the phase diagrams of solid 
hydrogen and its isotopes. We found that the Til phase 
line is interspersed by another phase, likely to be orien- 



tationally frustrated, as suggested by the experiments of 
Goncharov et al. (Ref . Hit) . For a thermal distribution of 
ortho-para rotors we find a reentrant phase diagram for 
both H2 and D2 . At low temperatures the order is short- 
ranged. Most experimental signatures of the long-range 
ordered phase are well known, 0, Q the short-range or- 
dered phase may be seen by probing the distribution of 
the local order parameter. While nuclear magnetic reso- 
nance .1?) is a useful probe, it may be difficult to apply 
at high pressure. The II' phase was found 14] by in- 
vestigating the Raman vibronic shift and it is known to 
be sensitive in ortho-para mixtures 0|. Rotational Ra- 
man lines are also sensitive to local order. Experimental 
studies to resolve the issues raised here would be helpful. 
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